Non-equilibrium phase transition in a periodically driven XY spin chain 
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We present a general formulation of Floquet states of periodically time-dependent open Markovian 
quasi-free fermionic many-body systems in terms of a discrete Lyapunov equation. Illustrating the 
technique, we analyze periodically kicked XY spin 1 /2 chain which is coupled to a pair of Lindblad 
reservoirs at its ends. A complex phase diagram is reported with re-entrant phases of long range 
and exponentially decaying spin-spin correlations as some of the system's parameters are varied. 
The structure of phase diagram is reproduced in terms of counting non-trivial stationary points of 
Floquet quasi-particle dispersion relation. 
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Introduction. Understanding and controlling dynamics 
of many-body quantum systems when they are open to 
the environment and driven far from equilibrium is an 
exciting and important topic of current research in the- 
oretical pQ [5] and experimental quantum physics |3] . In 
particular, since it has been recently realized that cer- 
tain emergent phenomena, such as quantum phase tran- 
sitions and long range order - previously known only in 
equilibrium zero-temperature quantum states [1] - can 
appear also in far from equilibrium steady states of quan- 
tum Liouville evolution j2l [5l-E]- In investigating dy- 
namical and critical many-body phenomena, quasi-free 
(quadratic) quantum systems play an important role as 
they are amenable to analytical treatment (see e.g. [H]), 
so many effects can be analyzed exactly or in great de- 
tail. For example, quantum phase transitions in non- 
equilibrium steady states have been observed either in 
quasi-free [SJ |S] , or strongly interacting [5] , or even dissi- 
pative [3 [TU] quantum systems in one dimension. 

At least two distinct types of quantum phase transi- 
tions with emergence of long range order have been iden- 
tified: (i) transitions which result from competition be- 
tween unitary and dissipative parts of the Liouvillian dy- 
namics [2J [7] , and (ii) transitions which occur in bound- 
ary driven systems |5l ^ (without any dissipation in the 
bulk) as a result some qualitative change in the proper- 
ties of the unitary dynamics (i.e. of the Hamiltonian) , 
say spontaneous symmetry breaking or bifurcations in 
quasi-particle dispersion relation i5^. Nevertheless, it has 
not yet been investigated whether such a structure of 
long range order can persist in the presence of external, 
periodic driving. Open and dissipative quantum system 
periodically driven by intense laser fields have been stud- 
ied by means of the so called Floquet theory [IT] , however 
most of the previous studies have dealt with essentially 
single-particle problems. From dynamical system's per- 
spective, it can be argued, that periodically driven or 
kicked quantum chains are probably the simplest and 
cleanest examples of quantum many-body ergodic dy- 
namics [12]. 



Furthermore, periodic time dependence offers new 
ways of encoding topological information (see e.g. [TBI). 
Periodically driven many-body systems have thus been 
proposed as candidate systems for studying topological 
order and detecting Majorana fermions in one dimension 
[13( I14j . On the other hand, quantum spin chains with 
particular specific Hamiltonians have been proposed as 
efficient ways for quantum information transfer |15) . En- 
gineering Hamiltonians and Liouvillians with particular 
properties in terms of a small set of primitive local opera- 
tions (realized by laser pulses) is one of the key problems 
in designing quantum computation, so it would be highly 
desirable to have simple and powerful techniques for the- 
oretical treatment of time-dependent open many-body 
systems, not necessarily time-periodic. 

In this Letter, we address the question of exact 
treatment of periodically driven quasi-free open many- 
body dynamics, which can be described either by time- 
dependent Lindblad equation, or equivalently, by discrete 
dynamical semigroups. We apply the idea of canoni- 
cal quantization in the Fock space of density operators 
[inillT] and solve for the covariance matrix of the (asymp- 
totic) stationary Floquet state in terms of the discrete 
Lyapunov equation. This yields an efficient general setup 
applicable to such class of problems in the language of 
classical control theory [TB]. Our method, which can be 
useful also in the case of non-periodic time-dependence, 
is then applied to treat an open Heisenberg XY spin 1 /2 
chain which is periodically kicked with a transverse mag- 
netic field. We find an appearance of two distinct phases 
of spin-spin correlations by changing the (bulk) system's 
parameters, namely the phase of exponentially decaying 
correlations and the phase of long range magnetic correla- 
tions. The phase diagram of the model is rather complex, 
but can be reproduced in terms of bifurcations of the Flo- 
quet quasi-energy quasi-particle dispersion relation. 

Stationary Floquet state. We consider a finite sys- 
tem of n (Majorana) fermions, which are described 
by 2?! anticommuting Hermitian operators Wj, j = 
l,...,2n, {wj,Wk} = '25 jk, which may be constructed 
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in terms of Pauli spin 1/2 operators o'^'^'^ via Jordan- with exphcit solutions 
Wigner transformation W2m-i = <^mY{m.'<m'^m'^'^2m = 
Y\m'<m ^m" ^6 shall be interested in a solution of a 
general time-dependent Markovian master equation for 
the system's density operator p{t), 



Q(f)=rexp - / dt'X(t') , 



(5) 



At 



p{t) ^ C{t)p{t), C{t)p=-i[H{t),p]+V{t)p, (1) 



where the dissipator may either be of a Lindblad form 
'^P = J2f,i'^L^pLl - {LlLf„p}), given in terms of (pos- 
sibly time dependent) Lindblad operators L^{t), or even 
of more general, Redfield form (see Ref. for the for- 
mulation in compatible notation). We focus on quasi- 
free dynamics where the Hamiltonian is given in terms 
of a quadratic form H = j, WjHj^kWk = w ■ Hw with 
antisymmetric imaginary matrix H and linear Lindblad 
operators — ^jlfi.jWj = If^ ■ w, G C^". Here 
and below, hatted symbols (such as the Liouvillian C) 
designate super-operators over 4" dimensional operator 
(Liouville) space, and bold-roman symbols (such as H) 
designate 2n x 2n or An x An matrices. 

Recently, a general approach to explicit analysis of 
dynamical properties of Markovian master equations of 
open many-body systems has been proposed [HI HEl HH] , 
which is based on quantization in the Fock space of (den- 
sity) operators, sometimes referred to as the 'third quan- 
tization'. In the context of quasi-free systems, if one 
is only interested in pair-correlations, dynamics of the 
covariance matrix Cj^k{t) ■— tT[wjWkp{t)] — 6j^k can be 
derived, either directly from ([T]) by observing canoni- 
cal anticommutation relations (CAR) among Wj^ or by 
means of fermionic super-operators over the Liouville 
space [51 [TB] 



C{t) = -X{t)C{t) - C(t)X'^(i) + iY{t). 



(2) 



X(t) = 4(iH(i) + Mr(i)) and Y{t) = 4(Mi(i) - MT{t)) 
are real 2n x 2?! time-dependent matrices, where M = 
Ml- + iMj is a bath matrix, given as M = l^ (g) l_^ 
for the Lindblad model, whereas the expressions for the 
Redfield model can be found in Refs. [HI [20] ■ In the 
time-independent case, one is interested in the steady 
state solution of (jij) C = 0, completely determined 
via the solution of the continuous Lyapunov equation 
XC + CX"^ = iY. For a general, time-dependent case, 
the general solution of ([2| for covariances can be sought 
for in terms of an ansatz 



C(t) = Q(i)C(0)Q^(0 - iP(i)Q^(t), (3) 



which results in two simpler equations, for 2n x 2n real 
matrix functions P(t), Q(i) 

Q(i) = -X(t)Q(i), Q(0) = 1, (4) 

P(t) = -X(<)P(<) - Y(i)Q-^(i), P(0) = 0, 



P(t) = -Q(t) / At'q,-\t')Y{t')q,-^{t'). 

Jo 

In this Letter we are interested in a time-periodic 
Markovian master equation ([l]), with £{t + t) = C{t), 
and looking for the stationary state pp, satisfying 

Zi(r,0)pF = PF, Uit",t') = texp dtCit)^ , (6) 

which we will refer to as stationary Floquet state (SFS). 
The covariancc matrix in SFS, Cp, can be given in terms 
of a solution of the discrete Lyapunov equation [18j 



Q(t)Cp-CfQ-^(t) =iP(T 



(7) 



which results from ^ after plugging C(0) = C(t) = Cp. 
It could be of some interest also to understand the struc- 
ture of the Floquet-Liouville spectrum of the many-body 
super-operator U{t, 0). We show in appendix to this Let- 
ter PTj that a complete spectrum of Z^(t, 0) can be writ- 
ten in terms of eigenvalues of the matrix Q(t). For ex- 
ample, the eigenvalue Ai of 2n x 2n matrix Q(t) of max- 
imal modulus (note that all eigenvalues of Q have mod- 
ulus less than 1), is also a maximal modulus eigenvalue 
of 4" X 4" matrix U{t, 0), and determines the time-scale 
t* — r/log |l/Ai| of relaxation to SFS. 

Kicked open XY spin chain. Let us now consider a 
special case of periodically kicked systems with time de- 
pendent structure matrix X(t) = X'^ + (5,-(i)X^ and time 
independent Y = SMj, where 6r{t) = r X^mgz "^''') 
is a periodic Dirac function with period r. Starting the 
one-period time interval just after the kick, we find ex- 
plicit solutions of Eqs. Q 



Q = exp(-rXi) cxp(-tX°), 
P = Q(exp(rX")Zexp(TX°)^ 



(8) 



Z), 



where 2n x 2n real antisymmetric matrix Z is a solution 
of the continuous Lyapunov equation 



X°Z + Z(X' 



0\T 



(9) 



Note that — iZ is the covariance matrix for the time- 
independent case TS^. In the following we apply our 
method to study the kicked XY Heisenberg spin 1 /2 chain 



Tl-l 



m = E 



1+7 
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-h6r{t) 



3 = i 



coupled to a pair of Lindblad baths at its ends, with four 



Lindblad operators Li/ 



^3A 



' ^ 1.2" n ■ 



Hence, the model is quadratic in Majorana fermions 
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FIG. 1. Residual correlations Crcs(T, h) (plotted in a log-scale, 
n = 64, with coupling constants Tf" = rf = 0.5, Ta = 0.3 
and T2 = 0.1) in SFS of kicked Heisenberg XY spin chain for 
anisotropy 7 = 0.1 (top left) and 7 — 0.9 (top right). Corre- 
lations that fall below 10~* threshold belong to white colored 
region and correspond to a phase with exponentially decaying 
correlations. At higher anisotropy LRMC phase ceases to be 
a simply connected region. Below we show some typical spin- 
spin correlation matrices chosen along the red line on 7 = 0.1 
phase diagram {n = 128). Leftmost panel clearly shows ex- 
ponential fall-off from the diagonal, the middle two belong to 
LRMC phase near the phase boundary, and the rightmost to 
a kind of antiresonance which occurs for h = p7r/(2r),p £ Z 
where the kick has no effect. 



Wj [!&, and the structure matrices can be written in 
2x2 block-matrix form as 

Yj- = 2irL ffyj,. + 2iT^ayS,.nSk 



(10) 



^ j .n^ k ,n 



where T 



L,R. -pLjR. 



± r^'^ and j, k — 1, . 



1. Us- 
ing standard linear algebra routines we have solved nu- 
merically the Eq. ([7| with (8[9) for different values of 
the system's parameters and different chain sizes up to 
n ~ 10'^. In particular we have focused on the exis- 
tence of long-range order in SFS of the model, so we 
have defined an order parameter - residual correlator 

Cr.s = (Ei^.T";="/'|C„.|)/(El^fc-"-"/'l). Inaccor- 
dance with the time-independent open XY model [5] we 
have discovered regions in parameter space (7, r, h) where 
Ci-es decays exponentially with n, and regions of parame- 
ters where we have correlations over large distances and 
Cros oc 1/n (the phase of long range magnetic corre- 
lations (LRMC)). However, in contrast with the time- 
independent case, we find here a very rich phase diagram 
(see Figjl]) with re-entrant LRMC phase. 

We find that qualitative properties of the phase dia- 
gram are completely independent of the bath parame- 




FIG. 2. Half the total number of non-trivial stationary points 
(as their number is always even) of quasiparticle dispersion 
0{it) for kicked XY chain (same parameters as in Fig.[l|. The 
structure evidently coincides with the phase diagram of resid- 
ual correlations Cics{T,h). 
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ters 2 , so we expect that theoretical explanation rests 
only on the bulk properties. In analogy to the time- 
independent case, where the phases have been explained 
based on quasi-particle dispersion relation of the infinite 
XY chain , we study dispersion relation for the Floquet 
modes Qu = e'^u in the absence of the baths (P = 0). 
Due to translational invariance of the infinite chain, the 
Floquet modes can be written in terms of the Bloch waves 
Uj = a^n) exp(ijK), namely T(K)a(K) = e'^(''^a(K), where 

T{k) = exp[-rXj,o] exp[-r(X?,oe"^ + X^.e'^'^)] (11) 
is a SU(2) matrix, with two quasi-energy bands 0i,2('t) = 



± arccos 



cos(2T/i)cos(2re(K)) 



sin(2T/i)sin(2re(As:)) 



cos K 



and e(K) — \J cos^ k + 7^ sin^ n is the quasi-particle en- 
ergy for the un-kicked XY chain, k e [— 7r,7r). As it has 
been argued in Ref. [S], LRMC phase is signaled by the 
appearance of a nontrivial stationary point k* in the 
quasi-particle dispersion relation, i.e. d0(K*)/dK = 
for K* ^ {0,7r}. In fact, near the phase boundary 
the nontrivial stationary points are emerging, say at 
K = 0, so the critical parameters can be obtained 
via Landau scenario by expanding the quasi-energy 
Q{k) = 6*0 + a(7, T, K)K^ + 6(7, t, h)K^ + ©(k^). Then, the 
nontrivial solution 7^ k* = ±-y/— 0(7, t, /i)/26(7, t, K) 
exists when a and h have different signs, thus a new 
pair of stationary points emerges when a{-^,T,h) = 
t(1 - 72) + i72sin(2r)sin(2r/i)/sin(2T(/i- 1)) = 0, 
which defines the critical field he by solving the equation 
a(7, T, /ic(7, r)) = 0. Such curves hc{j,T) exactly 
reproduce the rib structure of the phase diagrams in 
Figs. |1|2| and - where separating the regions between 
and 2 nontrivial stationary points - define the phase 
boundary. Furthermore, in Fig. |2] we observe that 
the phase diagram structure can be nicely reproduced 
by plotting the number of non-trivial stationary 
points K* of 0{k), the total number of stationary points 
being nj -f 2. This has interesting consequences on the 



4 




-8 -7 -6 -5 0.05 0.1 0.15 0.2 0.25 0.3 

In(Ah) Ah 



FIG. 3. Critical behaviour of correlations in non-LRMC phase 
(left) and LRMC phase (right) for kicked XY chain of size n = 
480, and r = 0.5, where A/i — h — he- In non-LRMC phase 
correlations decay as CrGs(?') oc exp(— r/^). We confirm, for 
two values of 7 = 0.1,0.5 that correlation length ^ diverges 
with critical exponent = 1/2 when approaching the critical 
point. In the right panel we show an abrupt transition of 
the order parameter Crcs when entering the LRMC phase. 
In the inset, zooming in the region near the phase transition 
point, we observe appearance of the correlation resonances 
(explained in the simplified context in Ref. ;20,)- Couplings 
Fj^'j are reduced by a factor of 10^^ w.r.t. Fig. 1 



temporal growth of block entropy in the operator space 
[22] which is conjectured as S{t) ^ |(n|j + 2) Int. 

One can define a critical exponent by approaching the 
phase boundary from the side where correlations decay 
exponentially Cj^k ^ exp{—\j — k\/£_), ^ ~ \h — hc\~''. Us- 
ing our heuristic quasi-particle picture, f can be derived 
from the imaginary stationary point wave-number near 
the critical point (7, r, /ic(7, ''')), k* ^ i\h — /ic(7, t)]"^^^, 
giving the critical exponent v — 1/2, which is confirmed 
by numerical simulations in kicked XY chain in Fig. [3] 
We have carefully checked numerically another signature 
of non-equilibrium quantum phase transition, namely 
the closing of the spectral gap A(n) = — ln|Ai|. We 
find clearly - similarly as in time-independent open XY 
model O [9] - that A(n) cx away from the phase 
boundary h ^ whereas A(n) oc for h — he- 
Conclusions. We presented a general setup for treating 
periodically driven open quasi-free many-body dynamics 
in terms of quantization in the Liouville-Fock space |23j . 
It is shown that covariance matrix for the stationary Flo- 
quet state satisfies a discrete Lyapunov equation, well 
known in classical theory of differential equations and 
in control theory [18] . Applying the method to study 
kicked open XY spin chain, we find a rich phase diagram 
of re-entrant phases of long-range correlation order and 
exponentially decaying correlations upon variation of a 
generic system's parameter. 

We acknowledge support by the grants Pl-0044, Jl- 
2208 of Slovenian research agency, as well Bessel research 
prize of A. v. Humboldt foundation. 
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SUPPLEMENTARY MATERIAL: 
FLOQUET-LIOUVILLIAN SPECTRUM OF 
DECAY MODES. 



Let us now briefly show how to construct the com- 
plete spectrum of decay modes of our quasi-free Floquet- 
Liouvillian theory. Following Refs. [HI [H] we de- 
fine An Hermitian Majorana maps a = {ai,a2), as 
ciij — {cj + Cj)/\/2, 0,2. j — i(cj — c^|)/-\/2, where canon- 
ical fermionic maps are uniquely defined by their action 

Cj\Pa) = 5ajfl\WjPa), c\\Pa) = 5a^s\WjPa) OVCr the 

complete Fock basis Pa — 2^"/^w"^ • • • W2n" '^^ ^hc op- 
erator space. Referring to Sect. 3.5 of 9 we find the 
Liouvillian propagator 



U{t,{)) = Ar(t)exp(ia- [lnU(t)]a) , (12) 



where, after a tedious but direct calculation, one finds 



"<')-'(Tc;5«,)s.s.^(;-), 

(13) 

with P(t),Q(t) given by (jsj and N{t) = [det Q(i)]i/2. 
Write shortly Q = Q(t), P = P(r), U = U(r), and 
assume that Q can be diagonalized 24 , Q = RDR^^, 



with D = diag{Ai, . . . , A2„}. Then, U is diagonalized as 
U = V-i(D®D-i)V, V= (^^q"' ^^'^f^S, (14) 

where V^^ = V-^J, J = a^®l2n- Now, in exact analogy 
with the time-independent case, we define normal master 
mode maps (6, & ) = Va, or 

6 = R"i(c+ Cpc^), h=Ii^c\ (15) 

satisfying CAR, {h^M = {b'^X} = 0, {h^.V^} = -^j^fe 
which diagonalize the many-body Liouville-Floquet map 

Z^(T,0)=exp |^|^(lnA,)S;.6,j . (16) 

As dynamics should be stable we have Vj, \ < 1. The 
complete set of 4" Floquet modes, U{t, 0) \]/) ~ A^ is 
then = Y\Y=i(^y^ \pf)i € {0, 1}, with decay rates 
A^ = Y\- AJ' . Note that the largest modulus eigenvalue 

Ai of Q determines the spectral gap of U{t,0), or the 
relaxation time to pp as t* — r/(— ln|Ai|). 

We note that the covariance flow ^ describes any, pos- 
sibly non-Gaussian, time-dependent state. Nevertheless, 
SFS Ipf) can be thought of as the right vacuum state of 
a non-Hermitian quadratic periodically time-dependent 
field theory, determined by bj \p-p) = 0, which implies 
that Ipf) is a Gaussian state, and all higher-order ob- 
servables can be expressed in terms of covariances Cp 
via the Wick theorem. 



